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Summary. — We present recent results on bulk observables and electromagnetic 
probes obtained using a hybrid approach based on the Ultrarelativistic Quantum 
Molecular Dynamics transport model with an intermediate hydrodynamic stage for 
the description of heavy-ion collisions at AGS, SPS and RHIC energies. After briefly 
reviewing the main results for particle multiplicities, elliptic flow, transverse momen- 
tum and rapidity spectra, we focus on photon and dilepton emission from hot and 
dense hadronic matter. 
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1. Introduction 



The investigation of nuclear matter under extreme conditions is one of the major 
research topic of nuclear and high energy physics. Experimental information about the 
properties of hot and dense strongly interacting systems is sought by analysing high 
energy collisions of heavy nuclei. To link specific experimental observables to the different 
manifestations and, eventually, phases of the strongly interacting matter, a detailed 
understanding of the dynamics of the heavy-ion reactions is essential. 

Numerous observables, such as hadronic and electromagnetic probes, their dynam- 
ical pattern and some specific correlations they exhibit, have been and are currently 
investigated in detail experimentally. All these observables are generally connected in 
a non-trivial way. A challenging task, in this respect, is a meaningful modelling of the 
heavy ion collision. Microscopic (transport) and macroscopic (hydrodynamical) models 
attempt to describe the full time evolution of the heavy-ion reactions and have played, 
in their various realizations, an important role in the interpretation of the experimental 
results over the last decades. 
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Recently, a third class of models, so-called "hybrid approaches" , has been developed. 
Hybrid approaches combine the advantages of transport approaches that are well suited 
to deal with the non-equilibrium initial and final states, with those of an intermediate 
hydrodynamic evolution, where, e.g., the equation of state (EoS) is an explicit input and 
phase transitions can be easily implemented. Such approaches were proposed 10 years ago 
[1, 2] and have since then been employed for a wide range of investigations [3, 4, 5, 6]. The 
hybrid approach discussed here is based on the integration of a hydrodynamic evolution 
into the Ultra-relativistic Quantum Molecular Dynamics (UrQMD) transport approach 
[7]. This integrated hybrid approach, called UrQMD v3.3 (^), has been applied by the 
Frankfurt group to the investigation of various hadronic observables in the broad energy 
range from Eiab = 2 — 160 GeV [8, 9, 10, 11, 12, 13] and, recently, of electromagnetic 
probes (photons and dileptons) [14, 15, 16, 17, 18], which have the unique feature of 
being sensible to the whole time evolution of the system. Note that in the hybrid model 
emission of virtual and real photons from the QGP phase can be explicitly accounted for 
provided that an EoS with a QGP phase is used for the hydrodynamical evolution. This 
constitutes a great advantage with respect to the pure hadronic transport models, where 
such emission cannot be easily implemented. In this proceeding, we briefly survey some 
of the results obtained both for the hadronic and electromagnetic sector. 

2. A hybrid approach to heavy-ion collisions 

Hybrid models generally schedule three dynamical stages: during the initial stage a 
microscopic transport scheme carries the incidentally colliding nuclei towards a stage that 
determines the initial conditions for the relativistic hydrodynamic equations of motion. 
At this stage one assumes local equilibration and the fireball is subsequently followed 
by a hydrodynamic evolution until the description is handed over to a final-state kinetic 
transport description. This late stage automatically furnishes a continuous freeze-out 
process, an important improvement compared to the otherwise employed prescription of 
an instantaneous freeze-out. Below, we want to briefly summarize the specific realization 
of the hybrid approach by the Frankfurt group and refer the reader to [7] for further 
details. 

During the first stage of the evolution the particles are described as a purely hadronic 
cascade within UrQMD. The coupling to the hydrodynamical evolution proceeds when 
the two Lorentz-contracted nuclei have passed through each other. At this time, the 
spectators continue to propagate in the cascade and all other hadrons are mapped to the 
hydrodynamic grid. Subsequently, a (3-1-1) ideal hydrodynamic evolution is performed 
using the SHASTA algorithm [19, 20]. The hydrodynamic evolution is later merged into 
the hadronic cascade. Two possible transition criteria procedures where tested. The first 
is the isochronous freeze-out (IF). In this approach, all hydrodynamic cells are mapped 
onto particles at the same time, once the energy drops below five times the ground state 
energy density in all cells. The second criterium is called gradual freeze-out (GF). In this 
approach transverse slices, of thickness 0.2 fm, are transformed to particles whenever in 
all cells of each individual slice the energy density drops below five times the ground state 
energy density. The employment of such gradual transition allows to obtain a rapidity 
independent transition temperature without artificial time dilatation effects [21]. When 
merging, the hydrodynamic fields are transformed to particle degrees of freedom via the 
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Cooper-Frye equation. The created particles proceed in their evohition in the hadronic 
cascade where final state interactions and decays of the particles occur within the UrQMD 
framework. 

3. Hadronic probes: a brief review 

Compared to treatments solely throughout by kinetic transport, as, e.g., by UrQMD, 
such hybrid strategies provide a well pronounced sensitivity of the transverse flow on 
the hydrodynamic part of the evolution [9, 11, 15]. Among others, they enhance the 
production of strange particles [8, 12]. In Fig. 1, the excitation functions of the A, S, J7, 
TT, K and p yields [22] are shown . The enhancement of the multiplicities for all particles 
with strange content compared to the non-hybrid results (dotted lines) can be clearly 
observed. 

Investigations of the longitudinal dynamics via rapidity distributions of various hadron 
species have shown that the latters are not too sensitive to the details of the dynamics 
for the hot and dense stage [8]. The rapidity distributions for, e.g., tt" and at three 
different energies (i?iab = 11,40 and 160A GcV) have been analysed in Ref. [8], where 
it was shown that the general shape of the distribution resulting in the hybrid approach 
is very similar to the one obtained by pure cascade calculations and in line with the 
experimental data [23, 24]. 

An additional outcome of such models is that elliptic flow is found to increase towards 
collider energies [13] compared to pure cascade calculations. 

4. Electromagnetic probes 

Electromagnetic probes, such as photons and Icpton pairs, are penetrating probes of 
the hot and dense matter. Once created these particles pass the collision zone essen- 
tially without further interaction and can therefore mediate valuable information on the 
electromagnetic response of the strongly interacting medium. 

In the hybrid approach, emission of real and virtual (dileptons) photons is treated 
as follows. During the locally equilibrated hydrodynamic stage the production of lep- 
ton pairs and direct photons is described by radiation rates for a strongly interacting 
medium in thermal equilibrium (so-called "thermal" photons and dileptons). In the evo- 
lution stage that precedes or follows the hydrodynamical phase, the emission is performed 
as typically done in cascade approaches. More specifically, dileptons emission is calcu- 
lated employing the time integration method that has long been applied in the transport 
description of dilepton emission (see e.g. [35]), whereas direct photon emission is cal- 
culated according to cross-sections for direct photon production from various channels. 
We will briefly discuss the main aspects of photons and dilepton emission in the hybrid 
approach here below. For further detail, the reader is referred to [14, 15, 18]. 

4'1. Photons. - The most important hadronic channels for the production of direct 
photons are tttt — >■ and irp — > [36], which both are implemented in the transport, 
as well as in the hydrodynamic phase. The cross-sections for cascade-calculations are 
taken from Kapusta et al. [36], while the rates used for the hydrodynamic description 
have been parametrized by Turbidc et al. [37]. Since no thermal partonic interactions 
are modelled in UrQMD, emission from a QGP-medium is only taken into account in 
the hydrodynamic part of the model. Several minor hadronic channels are only imple- 
mented in one of the two models, such as strange channels {e.g. Kn — > 7A'*) which 
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Fig. 1. - (Color Online) Excitation functions of the A, S, il, tt, K andp multiplicities (47r) in cen- 
tral (6 < 3.4 fm) Au+Au/Pb+Pb collisions. The results for the hybrid model with isochronous 
freeze-out (full lines), the hybrid model with gradual freeze-out (squares) and pure UrQMD- 
2.3 (dotted lines) are compared to experimental data (full symbols) from various experiments 
[25, 26, 27, 28, 24, 29, 30, 31, 32, 33, 34]. 



are only present in the hydrodynamic calculations, and ry-channels {e.g. irr] — > jn) which 
are only present in the transport calculations. Earlier investigations with this model 
have shown those channels to provide only minor contributions to the overall spectrum 
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Fig. 2. - (Color Online) Comparison of direct photon spectra (open squares) for central (0-20 %) 
and mid-central (20-40 %) Au-Au collisions at ^ysNN = 200 GeV [39] to different calculations: a) 
cascade (red crosses), b) hybrid model with HG-EoS (blue solid line), c) x-EoS (orange dashed 
line), and d) BM-EoS (violet dotted line). The contribution from initial pQCD-scatterings [40, 
39] has been added to all spectra. The spectra from central collisions have been scaled by a 
factor of 10^ to enhance readability. 



of direct photons. In the Quark-Gluon-Plasma, the rate used is taken from Ref. [38], 
where convenient parametrizations for the contribution of 2 -o- 2, bremsstrahlung- and 
annihilation-processes are given. The complete list of channels and a detailed explanation 
of the calculation procedure is provided in [14]. 

In Fig. 2 a comparison between direct photon spectra from hybrid model calculations 
and data from the PHENIX collaboration [39] for central (0-20%) and mid-central (20- 
40%) Au-|-Au collisions at ^/snn ~ 200 GeV is shown. The experimental data were 
obtained by extrapolating the dilepton yield to zero invariant mass [39]. The hybrid 
model calculations were performed using a hadron gas EoS describing a non-interacting 
gas of free hadrons [41] (HG-EoS, blue solid lines), a chiral EoS that follows from coupling 
a chiral hadronic SU(3) Lagrangian with a PNJL-type quark-gluon description [42] (%- 
EoS, orange dashed lines), and bag model EoS that exhibits a strong first order phase 
transition between a Walecka type hadron gas and masslcss quarks and gluons [20] (BM- 
EoS, violet dotted lines). Pure cascade calculations are indicated by red crosses. All 
calculated spectra include the (A'^coii)-scaled prompt photon contribution taken from [40, 
39]. We observe that in both centrality-bins, the direct photon spectra obtained with the 
BM-EoS and x-EoS, which include a phase transition to a deconfined state of matter, 
are significantly higher than the hadronic HG-EoS-calculations. Similar enhancement of 
photon production due to QGP emission was found already at SPS energies [14]. 

4'2. Dileptons. - Invoking vector meson dominance the emission rate of Icpton pairs 
from a strongly interacting medium can be related, at low invariant masses, to the spectral 
properties of the vector mesons, with the p meson giving the dominant contribution. The 
thermal dilepton rate reads then [43] : 

(1) = -^^^fBilo; T) Im D,{M, q; T, ^b) , 
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where a denotes the fine structure constant, = Qq ~ 1^ the dilepton invariant mass 
squared, /b the Bose distribution function (for a moving fluid this must be substituted 
with the Jiittner function), and L{M'^) a lepton phase space factor that quickly ap- 
proaches one above the lepton pair threshold. The electromagnetic response of the 
strongly interacting medium is then encasted in lmDp(M,q]T, hb), the imaginary part 
of the in-medium p meson propagator, 

(2) Dp{M,q-T,iiB)^jT-^ 2 i (M 

-m^p- T.p{M, q; T, (xb) 

Over the years, strong evidence has been accumulated pointing to the conclusion that 
the inclusion of the in-medium contributions to the p meson self energy is mandatory 
for a proper description of the low invariant mass region of the dilepton spectra (sec e.g. 
Ref. [44]), typically characterized by the emerging of an enhancement with respect to 
the standard hadronic cocktail [45]. A nice aspect of the hybrid approach is that the 
hydrodynamic stage allows for a transparent inclusion of the in- medium spectral function 
of the vector meson, conceptually problematic in transport calculations, in analogy to 
fireball model calculations [46, 47]. An advantage with respect to the latters is the use 
of a dynamical model for the description of the heavy-ion collisions, which might help to 
shed some light on dynamical aspects as dilepton radial flow or similar. 

In this application, the self-energy contributions taken into account are Sp = + 
Y,pn_^_j^pN^ where T,° is the vacuum self-energy and E'"' and denote the contribution 
to the self-energy due to the direct interactions of the p with, respectively, pions and 
nucleons of the surrounding heat bath. The self-energies have been calculated according 
to Ref. [48] , where they were evaluated in terms of empirical scattering amplitudes from 
resonance dominance at low energies and Regge-type behaviour at high energy. 

In Fig. 3 hybrid model calculations are compared to recent acceptance-corrected NA60 
data [49] . The calculations have been exemplary performed using a hadron gas equation 
of state for the hydrodynamical evolution. This is enough to point out qualitatively the 
main features of the present approach. More systematic studies and discussions will be 
presented elsewhere [50]. We observe that the cascade emission dominates the invariant 
mass region around the vector meson peak for both low and intermediate transverse pair 
momenta pt- At low pT (left panel of Fig. 3), the very low invariant masses, M < 0.5 
GeV, are flUed by the thermal radiation with in-medium spectral function. The sum 
of both contributions, however, leads to an overestimation of the vector meson peak 
region at low transverse momenta of the dilepton pair. The reason for the discrepancy 
might partially lie on the speciflc spectral function used here and/or, presumably more 
severely, on the eventual presence of not yet negligible residual in-medium modiflcation 
of the p meson spectral function during the cascade stage, that are here neglected. With 
increasing pT (right panel of Fig. 3), dilepton emission at very low invariant masses is 
reduced and the total spectra are almost completely determined by the cascade emission. 
Indeed, thermal emission had already shown discrepancies for px > 1 GeV [46] , pointing 
to the necessity to account for non-thermal contributions. In the present approach, the 
latters appear quite naturally. 

Finally, we would like to mention that, in analogy to the direct photons calculations, 
the present approach easily allows for the inclusion of dilepton emission from QGP, which 
is thought to play a role in the intermediate mass region of the dilepton spectra. 



* * * 



HADRONIC AND ELECTROMAGNETIC PROBES OF HOT AND DENSE MATTER IN A BOLTZMANN+HYDRODYNAMICS MODEL OF RI 



p.|-<0.2 GeV 

HG EOS 



I I .1 . L_L 



• NA60 

— - hydro emission 

— cascade emission 

— sum ( p) 




3.2 0.4 0.6 0.8 1 

M [GeV] 



1.6<p.|-<1.8 GeV 



"T" 



"T" 



• NA60 

— - tlydro emission 

— cascade emission 

— sum ( p) 




Fig. 3. - (Color Online) Left panel: Acceptance-corrected invariant mass spectra of the excess 
dimuons in In-In collisions at 158A GeV for transverse pair momenta pr < 0.2 GeV, compared 
to hybrid model calculations based on thermal radiation from in-medium modified p meson 
spectral function (dotted-dashed line) and non-thermal cascade emission (dashed line). The 
sum of the two contributions is depicted by the full line. Experimental data from Ref. [49]. 
Right panel: Same as in the left panel, but for the transverse momenta window 1.6 < px < 1.8 
GeV. 
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